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Abstract. This article presents a squeezing transformation for quantum 
systems associated to finite vector spaces. The physical idea of squeezing here is 
taken from the action of the usual squeezing operator over wave functions defined 
on a real line, that is, a transformation capable to diminish (or enhance) the mean 
square deviation of a centered distribution. As it is discussed, the definition of 
such an operator on finite dimensional vector spaces is not a trivial matter, but, 
on the other hand, has obvious connections with problems such as spin squeezing 
and (finite) quantum state reproduction. 



PACS numbers: 03.65.-w,03.65.Ca,42.50.-p,42.50.Dv 

1. Introduction 



Quantum mechanics has a lot of peculiar features. One of them is that it is able to cope 
with physical quantities which manifest themselves in two quite different ways, namely 
those which are better described by continuous variables (a case which we hereafter 
shall refer to as the 'continuum'), and those described by finite and discrete variables 
(mentioned as the 'discrete'). These both 'versions' of quantum mechanics have their 
own idiosyncratic features, and, despite the late rise of interest in finite dimensional 
systems triggered mostly by the field of quantum information (which in fact deals 
with both continuous and discrete cases), one can safely say that the continuous case 
have a much richer (and better established) variety of theoretical tools. Among the 
various examples one could name here, a particularly sound one is that of squeezing. 
The importance of the squeezing transformation (and the squeezed states) in general 
quantum mechanics, quantum optics and, obviously, in quantum information theory 
associated to continuous variables goes literally without saying, and references on the 
subject are countless (some good examples, among many, of textbooks and review 
articles are [1] [2]). On the other hand, however, even the idea of squeezing is hardly 
touched upon in regard of finite dimensional quantum systems. 

In the continuum, (that is, in the quantum mechanical description of the 
usual, continuously labeled, canonically conjugated variables of position and linear 
momentum) the unitary transformation of squeezing acts reducing the dispersion 
related to, for example, the position observable (for zero- mean distributions). 
The Fourier transform connection between the wave-function in the position and 
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momentum representations prevents from any conflict with the conceptual content of 
the indeterminacy principle, as the same transformation 'stretches' the wave function 
in the dual momentum representation. The domain of such wave-functions is dense 
and infinite, and this both attributes are crucial in order to be possible to squeeze any 
given wave-function. In rough terms, if the domain were not infinite, from where the 
points (in the new, squeezed, representation) near the boundaries of a given region 
(that can always be made arbitrarily large) would have come from? And if the domain 
were not dense, to where all this points would go? In other words, how can you squeeze 
a given wave function, if this wave function is defined over a discrete and finite grid, 
keeping the grid fixed? All those "ifs" are obviously unnecessary concerns as long 
as one is dealing only with the usual Cartesian coordinates, but lie in the core of 
the problem if one would like to consider the problem of squeezing in discrete finite 
dimensional quantum systems. 

A well known squeezing operator in finite dimensional quantum systems would 
have immediate consequences in the quantum information theory for such systems. 
The no-cloning theorem [3] asserts that quantum states may only be completely 
reproduced if and only if they are an element of a given known vector basis. Naturally, 
the higher the overlap of a given state with some known basis element, the better 
such state can be cloned [4]. Now, as we will discuss, a squeezing operator can be 
regarded as an operator which reduces (or increases) the dispersion (the mean square 
deviation from the mean) of any given state with respect to a given observable. In 
other words, such operator increases (reduces) the overlap between a given state and 
an untransformcd clement of a basis set, thus enhancing (dimishing) the fidelity of 
the reproduction of such state. Naturally, the actual implementation of any unitary 
transformation in a real device have its own difficulties and limitations, which always 
have to be taken into account. However, it is always better to face a technical difficulty 
(the greater it may be) than a theoretical impossibility. 

The aim of this article is thus to develop a mathematically consistent approach 
for squeezing in finite dimensional quantum systems, introducing what might be 
called squeezed discrete oscillator states and, naturally, the squeezing operator. The 
unitarity of the squeezing operator here proposed is a delicate issue. To put it in 
a paragraph, first a provisional squeezing operator is presented, which is a genuine 
squeezing operator in the sense that when properly applied it may reduce (or enhance) 
the mean square deviation of any finite distribution. Such operator, however, is not 
unitary. Unitarity is reached only in approximate fashion, so an objective criterion of 
approximate unitarity is introduced in an appendix. Therefore, unless explicitly stated 
otherwise, all use of the term "unitary" throughout this paper is in the approximate 
sense there defined. It must be stressed that the approximations presented arc 
extremely good. 

In the forthcoming section II the difficulties above presented regarding squeezing 
for finite dimensional quantum systems are discussed, and necessary initial discussions 
in order to achieve a finite dimensional squeezing operator are made, including a brief 
outlook of known properties of squeezing in continuous systems, which is absolutely 
necessary to guide the discussion. Section III is devoted to the (wayward) definition 
of the squeezing operator. The actual squeezing is presented in section IV, while final 
remarks and discussions are left to section V. 
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2. The theoretical problem of squeezing in finite dimensional quantum 

systems 

In quantum mechanics of usual, infinite dimensional and continuously labeled Hilbert 
spaces, the squeezing operator acts over the position eigenstates as 

S(X)\x) = -±=\kc), (1) 

where A is a real number which parametrizes the transformation. It is not hard to 
see the scaling effect that this transformation has on a given wave function, as a state 
represented by ip{x) would be transformed in to one represented by ip(\~ 1 x) 1 thus 
the name squeezing (stretching, if A is larger than one). However, in a finite, discrete 
domain, an operator which maps a state |fc) onto |Afc) (with k some discrete label of 
a given finite family of eigenstates) hardly could be obtained or even interpreted as 
a squeezing operator. In such domain, if A is larger than one, Xk would undoubtedly 
fall off the domain of labels at least for some values of k, and if A is smaller than 
one, as the domain is not dense, Afc would not be in the domain for some values (even 
all) of k either. A more sophisticated possibility would be to consider the product 
\k as a modulo N operation, as considered in [5]. In this case the difficulties just 
pointed out surely do not arise (at least for integer A and prime N). However, such 
a transformation cannot be interpreted as squeezing. It is a well known fact from 
number theory that, given a set of integers D obeying {k € D\0 < k < N — 1} (which 
is a complete set of residues modulo N), the set {Afc (modiV)} , for integer A, coincides 
with D, but in some different order [6]. In other words, this is rather a shuffling 
transformation than a squeezing (or stretching) one. The shape of an arbitrary wave 
function would be complete and radically altered by such transformation, and, more 
importantly, the effect of this transformation on the dispersion would depend heavily 
on the initial wave function. 

Thus the action of the squeezing operator over the families of eigenstates of 
position and momentum gives no valuable hint of how an equivalent operator in a 
discrete domain would act, and we must look elsewhere if we want to find a good 
starting point for the discrete version of squeezing. 

2.1. Oscillator states, width and Fourier transforms in the continuous 

The eigenstates of position and momentum are not the only way to look at the problem 
of squeezing. In a Schroedinger inspired approach to the issue, in this sub-section 
we look at some simple (and well, well known) properties of the harmonic oscillator 
states. This seemingly useless (but brief) repetition of general knowledge will prove 
itself more than pertinent in the next sections, as all expressions there obtained will 
end up being discrete analogs of the ones here presented. 

The ground state of the harmonic oscillator is described by a Gaussian function 
of unit width (in natural units), 

X), W° = TT 1 ^ 

written in terms of the position eigenstates \x). The whole set of states is defined 
through the Hcrmitc polynomials, 

\x), M n = V¥n~t.Tr 1/i . 



I f 

|0 > = W_, 



dxexp 



dxH n (x) exp 
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In order to account for different widths of the harmonic oscillator wave functions, it 
is possible to define 



|n;A) 



1 



N n y/XJ- 



dxH, 



(5 



exp 



I 

2A2 : 



\x). 



In the momentum representation, the same states are represented by the (integral) 
Fourier transforms of the functions above 



|n;A) 



dpH n (pX) exp 



A 2 



-P 



\P), 



(2) 



which have their widths values inverted in the process. 

In an algebraic sense, if the Fourier operator maps a position cigcnstatc on a 
momentum eigenstate as 

T\x) = \p), 

it directly follows that 

F\n;X) =i n \n;l/X), (3) 

and the oscillator states with unit width are eigenstates of J 7 . By its turn, the squeezing 
operator acts over the unit width ground state as 



S(A)|0;1) = 



dx exp 



1 

2A2 : 



\x) = |0;A), 



aa°Va. 

altering the width for the value A. In fact, in general 
S(A)|n) = |n;A). 

and, as the parameter A controls the variance associated to all oscillator states, the 
state is therefore "squeezed" (for A < 1). One must not forget that such "squeezing" 
may actually be a stretching if one is looking at the momentum representation, as 
shown by Eq.(2). 

Finally, oscillator states are complete and orthogonal, 



(n'; A|n; A) = 5 n , n * ^ \n; A) (n; A| = i, 

n=0 

what allows us to write the squeezing operator as 

00 

S(A) = S(A)5>;l)(n;l| 



n=0 



S(A)= X>;A)<n;l|. 



(4) 

(5) 
(6) 



n=0 



2.2. Discrete oscillator states 

2.2.1. Preliminaries Let the set {\uk)Y k= _ t be some complete and physically 
meaningful set of states [7, 8], spanning a N dimensional vector space QJ, associated 
with a given physical quantity O, so that 

0\u k ) = o(k)\uk), 

where O is the Hermitian operator representing the physical quantity, with its 
associated set {o(fc)}^, = _ £ of eigenvalues standing for the possible outcomes of physical 
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measures of O. Denning a "coordinate" operator for such discrete set of states as 
Q|ufe) = k\u k ), the dispersion a associated to any state may be written as 

a=(Q 2 )-(Q) 2 . 

If a given physical system is in a state \uj Q ) associated to the value o(jo) of O, the 
average value of the "coordinate operator" Q is obviously jo and the dispersion is 
zero, that is, the distribution is sharply peaked at jo- 
Let us also introduce, with the aid of the Discrete Fourier Transform (DFT), the 
complementary set 

i 

\vk) = exp 

such that a given state \ip), represented by a ip(k) wave function in the {|«fe)} 
representation, is represented by the DFT 4>{k) in the {|wfc)} representation. It is 
not hard to see that these two representations play in the discrete the role played by 
{|x)} and {\p}} in the continuum [8]. 



2iri 
'~N~ 



kj 



2.2.2. Width parametrized eigenf unctions of the DFT Now we leap in to much lesser 
known grounds. Reference [9] presents some very eloquent results regarding the 
discrete Fourier transform (DFT) of Jacobi Theta functions, that completely parallel 
the ones presented in subsection II. a. In particular, there it is shown that the set of 
functions 



fn(j,0 

behaves, under DFT, as 

JV-l 



' N d r ' 



t=o 



Vn 



3=0 



2ni ' 



fn(j,Z)=i n fn(k,C 1 )- 



The Jacobi ^3-function explicit form is [10] 

00 

#3 (z, t) = ^ exp [i7rTa 2 ] exp [2niaz] 



(7) 



(8) 



and with the help of the Hermite polynomials generating function, the f n (j, 
functions may also be put in the equivalent forms 



/n(j;0 = exp 



H n ( -(aN + j) 



fn{j;0 = 



00 

■i) n Y ex P 



a — — 00 



ir , 2iri 
—£ 2 a 2 H ja 



These functions are periodic, that is, f n (j + N; £) = f n (j; ^). As the second argument 
of the Theta functions controls their width [9], expression (7) ensures that the 
parameter £ controls the dispersion associated to those functions, as it is also evident 
from figures 1 and 2. The DFT preserves the analytic form of the functions / n (j;^), 
mapping £ into that is, a given width value into its inverse, in close analogy to the 
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behavior of the Gaussians under integral Fourier transforms. The f n {j',0 functions 
are symmetric about zero, and thus have zero mean in the sense that always 

e 

£i l/n(j,£)l 2 = o, e = ^. (9) 

Therefore, one may define discrete oscillator states as 
1 1 

\ n ^) = T^ E /«( fc ;0K>. ne[0,JV-l], (10) 

where A/"„^ is the appropriate normalization factor. From equation (9) it follows that 
(n;£|Q|n;£) = 0. By stating that "the parameter £ controls the dispersion" we mean 
that a, when calculated with respect to a discrete oscillator state, is a monotonically 
increasing function of £. It is evident from the definitions that the "coordinate" 
operator Q is NOT diagonal in the discrete oscillator states representation. 

Although it is obvious, it must be kept clear that the only analogy used here 
is the one between the discrete oscillator states and the usual continuous harmonic 
oscillator eigenstates, and NOT general number states in a Fock space. 

Algebraically, if the Fourier operator in the discrete domain is Td\uk) — \ v k), it 
is then immediate to verify that 

^|n;0 = Kr 1 ), (11) 
in analogy to Eq.(3). And, equivalently, 

1 1 

K0 = T7 E fnfcC^Vk), ne[o,N-i\, 

that is, the state |n;£) is represented by / n (fc;£ _1 ) m t nc dual {|«fe)} representation, 
again in a perfect analogy with the continuous case. 

It might be safely assumed that, for odd N, the first N discrete oscillator states 
are complete[ll]. For even N, the / n (fc; 1) functions have some interesting (and still 
not fully comprehended) properties, and the function associated to the upper bound 
N — 1, /jv_i(fc; 1), is identical to /jv-s(fc; 1) up to a minus sign. In view of that, for 
even dimensions we must substitute /jv-i(fc;l) by /jv+3(fc; 1) (as i N ~ x = i N+3 ) in 
order to have a complete set. In the following, for simplicity, we restrict our notation 
to the odd case. 

A little digression on the unit width states {|n; 1)} is necessary. The first 
important question to address is if this set is complete. First, in reference [11], Mehta 
conjectures that this set of states is linearly independent (and therefore complete). 
Such conjecture is heavily supported by numerical evidence f. We present the actual 
values of the overlap matrix (n 1 ; l|n; 1), for the case N = 13, which arc well illustrative. 
As those states are eigenstates of the Fourier operator with eigenvalue i n , (n; l|n'; 1) 
is exactly zero if i n ^ i n . Non-zero off-diagonal matrix elements occur for \n — n'\ 
a multiple of 4. In fact, in general, only for the larger values of the indices n and n' 
the values of the off-diagonal elements become significant [12]. Reference [9] presents 
a detailed study of the inner product (n; l\n'; 1). If the width £ is different than 
1, there arc more non-zero off diagonal elements. However, for values of £ in the 

X In reference [9] there is a misjudgmcnt of Mchta's conjecture. Guided by his own concerns on 
the subject, the author assumed that Mehta has conjectured on the orthogonality of the / n (A;;l) 
functions, while the actual conjecture is about the linear independence of the set. 
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Figure 1. Plot of the functions fo(k;l) and /o(fc;1.3). As those functions are 
defined on integers, the dotted lines only guide the eye. Two characteristics are 
explicit: First, the resemblance with the continuous harmonic oscillator functions, 
and the increase in width when one goes from the squeezing parameter 1 to 1.3. 
The parameters are choosen so that the squeezing effect is visible. 



vicinity of 1 (0.9 < £ < 1.1) the deviation from the values obtained for £ — 1 is 
small. In fact, for the lower values of the index n, the inner product (n;£|n';£) does 
not (significantly) depend on £, and therefore a transformation mapping \n'] 1) on to 
|n';£) is (approximate to) an unitary one, a fact which is duly explored in [13]. A 
precise notion of "approximate unitarity" is presented in the appendix. 

3. The squeezing operator for finite dimensional quantum systems 



The results of the previous section show the deep analogy between discrete and 
continuous oscillator states, as both families are parametrized by its width value, which 
is inverted under Fourier transformation, not to mention the shape of the functions 
itself, as can be seen in figures 1 and 2. In fact, one may obtain the discrete / ra (fc;£) 
from the continuous oscillator functions through the Poisson summation technique 
[14]. Such results give us a path to follow in order to define a squeezing operator 
suitable for finite dimensional quantum systems. 

3.1. Preliminary attempt 

The general idea is that the eigenstates of the harmonic oscillator provide a "natural" 
representation of the squeezing operator, and that discrete oscillator states provide 
a discrete analog of that representation on finite dimensional quantum systems. So, 
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Figure 2. Plot of the functions /i(fc;l) and /i(fc;1.3). The dotted lines only 
guide the eye. 



bearing in mind expression 6, we can define an finite dimensional squeezing operator 
of the form 

N-l 

S p (O = X>;0<n;l|. (12) 

n=0 

where the subscript p stands for 'provisional'. Now, let us suppose for a moment that 
the set {|n; £)} is orthogonal for all values £, or at least for a reasonable range around 
1. Such an operator, in that case, would perfectly parallel the continuous squeezing 
operator. By construction, it would squeez (stretch) a |m; 1) state characterized by 
width 1 to a smaller (higher) value £, that is 

S p (0|n;l> = |n;0, 

just like in the continuous case. Also, it would follow (after some algebra) that 
"p(0 — 3 p (£ _1 ), what, if Hp is unitary, leads to the convenient property 

s p (r 1 ) = s P - 1 (o- 

The unitarity of the 3 P (£) operator would follow directly 

N-l N-l 

V V |n;0 (n;l|m;l) (m;£| 

n=0 m=0 

O m , n 

N-l 

J2 \n;t)(n;t\ = t, 

n=0 

N-l N-l 



sJ(0s P (0 = 
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Table 1. Overlap matrix for N = 13 and 5 = 1- Values presented with two 
decimal digits were numerically calculated. The blocks were chosen so as the 
first diagonal block is the identity matrix, and the off-diagonal blocks have only 
elements smaller than 0.01 



N-l 

St(0E p (0= £ M)<n;l| = i, 

where the completeness of the {|n;£)} states is explicitly used. 

Under the orthogonality assumption, then, it is true that the above operator 
is unitary, and squeezes, or stretches, the discrete oscillator states, which do form 
a complete set in which any other state might be expanded. The squeezing, or the 
stretching, is understood in the sense that the dispersion a is diminished, or enhanced. 
However, as stated above, the discrete oscillator states are, in general, not orthogonal. 
Therefore, we must amend our initial (and naive) approach. 

3.2. First attempt: Non-orthogonal decomposition 

First of all, as the set {\n;£,)}^~Q is not an orthogonal set, in order to work properly 
with it we must adopt the usual techniques valid for non-orthogonal basis sets [15]. 
Due to the non-orthogonality, the operator defined on 12 does not map a unit width 
discrete oscillator state in to another discrete oscillator state of different width. This 
must be accomplished by means of a dual set with the property 

(m;£|n;£) = <5 m , n . 

In this form the closure relation may assume one of the forms 

JV-l N-l 

5>;0(»;£l = £l n ;0 M = i. 

so the decomposition of a given state would read, for example, 

N-l N-l 

IV) = £ |n;£)MV>> = £v4 \n;0- 

n=0 n=0 
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It is worth emphasizing that such a dual set always exists if the overlap matrix defined 
by the elements (j; £\l; £) has an inverse. Armed with the above results one could study 
the properties of the given operator 

s(o = X>';0(n';i|- 

n'=0 

Now, the action of this operator over the discrete oscillator states is clearly E(£)|n; 1) = 
|n;£). Therefore, in fact it squeezes, or stretches, the discrete oscillator states. It is 
interesting to note that the operator 

N-l 
n'=0 

satisfies the relation 

3(03(0 = s(0s(0 = i, 
thnss(0 = s- 1 (0- 

The condition for unitarity would then be S^(0 = 3(0, which in general does 
not hold. 

3.3. Second attempt: Squeezing conjecture and approximate unitarity 

The non-orthogonality of the states is a serious hindrance for the construction 

of an unitary squeezing operator, but it can be at least approximately circumvented. 
From now on, we follow the criteria defined on the appendix for "approximate 
unitarity" , and all further references to "unitarity" are made in this approximate 
sense. 

Some basic definitions are necessary: For a given dimensionality N, let Vf 
be the subspace spanned by the set of the first Ni (where Ni shall be defined 
a posteriori) discrete oscillator states, {|m; Olm^o^" 1 > an( ^ tne "remaining" 
subspace (orthogonal to V?) so that V, Vf = V. Accordingly, 1; and lh are the 
respective identity operators on those subspaces. 

In that case we can then define an unitary operator of the form 

JV,-1 

s„(0= E |n;0(«;i| + ifc. (13) 

n=0 

where i/j = X)n=A l n ' CI- Now, such an operator still has a lot of features which 
parallel those of the continuous squeezing operator. By construction, it squeezes 
(stretches) a |ra;£) state characterized by width 1 to a smaller (higher) value £, that 
is 

s„(0l«; i) = 

provided that n < Ni (otherwise it leaves the state unchanged). It is not hard to see 
that £„(£) will be unitary as long as 

(m; l|n; £} ~ 0; for m < Ni < n. 

In order to put forth this approximation, for a given dimensionality N of the 
original state space, one identifies, as in done in table I, two diagonal blocks of the 
overlap matrix (n;£|n';£), where the lower block is the identity matrix. The blocks 



A squeezing formalism for finite dimensional quantum systems 



11 



must be chosen so that all the elements in the non-diagonal blocks are negligible. The 
first block dimensions are Ni x N[, whereas the remaining block is x N^, obviously 
with Ni + N h = N. The appendix gives more details on this issue, but the general 
idea is the one described above. 

We are now in position to make a squeezing conjecture: The unitary operator 
which squeezes the discrete oscillator states belonging to vf squeezes all states defined 
on that subspace. 

4. Actual Squeezing 

Let us then briefly see the action of the finite dimensional squeezing operator. As the 
set {\n;{;)}n~Q is complete, an arbitrary state can always be described by 

N-l 

k>> = 5>(»;0|n;0- (w) 

One can, in principle, arbitrarily choose the value of £ in decomposition (14). Choosing 
then £ = 1, the action of H u (£) over \ip) may be simply written as 

N-l N,-l N-l 

S(0\<p) = 3(£) Y, <P{n; l)|n; 1> = £ <P{n; + E ^ ^ X >- 

n=0 n=0 n=Ni 

In figure 3 we present the actual squeezing of a square wave function. A number of 
different examples could have been shown, but their behaviour is basically the same. 
Once the distribution is roughly centered around zero, and can be reasonably described 
by the lower discrete oscillator states, the effect of the squeezing transformation on 
it will be appreciable. In the particular case here shown, the squeezing effect can be 
easily grasped from the figure. 

In fact, although the squeezing formalism here presented is clearly built to 
work for zero-mean distributions, one could easily generalize it simply by adopting 
conveniently shifted discrete oscillator functions, or making use of the well known 
unitary displacement operators for finite dimensional systems [7, 16]. 

The effectiveness of the result shown here should not be underestimated due to 
the presentation of only one example. Any centered, bell-shaped distrution would 
be even easyer handed by the discrete squeezing operator, the square wave function 
was choosen exactly because it is very different from the discrete oscillator states 
(and the discrete squeezing operator, by definition, correctly squeezes those states). 
The squeezing operator is able to unitarily control the dispersion of any centered 
distribution. 

5. Final remarks 

A straightforward analogy with the continuous does not give always immediately the 
best possible direction to follow in the discrete, as sometimes the peculiarities of such 
scenario should be explicitly taken into account. The continuous is a limiting, and 
thus particular, case of the discrete [8, 17, 18], and not necessarily all the conveniences 
of the continuous formalism are to be found in the discrete approach. However, such 
conveniences (that is, established theoretical tools with a commonly accepted physical 
interpretation) are often desirable, and a lot of effort is usually made to achieve them 
in the discrete finite dimensional scenario (a good example is the discussion about 
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Figure 3. Square wave function before and after squeezing transformations of 
parameters 1.1 and 0.9 



discrete Wigner functions). As we argued above, in the search for a finite dimensional 
squeezing operator, the position-momentum representation of the squeezing operator 
is misleading, and another starting point had to be found. 

As we have seen, such starting point is a complete set of states which is 
parametrized by its width. The idea can be summarized in the following lines: We 
construct an operator which maps unit width discrete oscillator states on to discrete 
oscillator states of smaller (higher) width. As discrete oscillator states are complete, 
we may conjecture that such operator maps any given state on to another of smaller 
(higher) width, thus characterizing squeezing (stretching). 

The first attempts to produce such operator, however, came up with something 
which is not unitary, as the discrete oscillator states lack orthogonality. We have seen, 
however, that if one works only with the lower discrete oscillator states, unitarity can 
be reached in approximate fashion (in other words, the squeezing transformation is 
approximately unitary in the sub-space of the lower discrete oscillator states). This 
means that, as states |0; 1) has been successfully used to define a reference state [13] 
for finite dimensional coherent states, due exactly to its convenient properties under 
DFT [19, 20, 21], the squeezing transformation over such reference state can be safely 
assumed to be unitary. 

In the above cited reference, in fact, the properties of finite coherent squeezed 
states, obtained as finite displacements of a squeezed reference state [0; £) are studied, 
and the first applications of the here presented approach in quantum information 
theory (which may well be the natural ground for applications of this formalism) 
are shown. Other interesting field in which this formalism could be immediately 
implemented is spin squeezing [22]. 

One could imagine that an alternative way of circumventing the non-unitarity of 
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the discrete oscillator states would be to obtain a new, orthogonal basis by means 
of some orthogonalization process, like the Gram-Schmitd procedure. Three such 
attempts have been extensively tested and none of them has proved to be satisfactory, 
all for the same reason. The new, orthogonal, sets of states, obtained trough 
all the attempted orthogonalization processes, no longer have their dispersion as 
a monotonically increasing function of the original squeezing parameter, so such a 
formalism, although mathematically consistent, looses it physical significance. One of 
the orthogonalization processes attempted was the Gram-Schmitd itself, the second 
was a variation of it, and the third one was a method which is based on the 
diagonalization of the overlap matrix. 

One must stress that key feature of squeezing, aside from its scaling, shape 
preserving effect on the wave function, is rather the consequence of this scaling, which 
is to diminish (or enhance) the mean square deviation of a given variable (in the case of 
a zero mean distribution). In fact, if taken together with a displacement operator, the 
squeezing operator may squeeze any distribution around its mean. The formalism here 
presented is focused precisely at this point (our squeezing operator could be rather 
called root mean square deviation diminishing operator, but, as Bohr once stated, 
precision is the complementary counterpart of clarity), and it ultimately gives one 
a procedure with which one may reduce the dispersion associated to any quantum 
variable defined on a discrete domain. 

As it is well known from quantum information theory, absolute fidelity of a 
quantum teleportation protocol, if one does not want to destroy the original system, 
can only be achieved if the original state is a member of an orthogonal basis set. 
For a general state, it is clear that, the higher the overlap of this state with a given 
basis element, the higher the fidelity of the process. The squeezing transformation is 
completely capable to enhance, or diminish, such fidelity. Therefore, one could argue 
that the ultimate end of the here proposed finite dimensional squeezing operator may 
be to squeeze out quantum information to the last q-bit. 

Acknowledgments: The author thanks K. Zan and J. Bclthcr for much needed 
help with numerical calculations, and M.A. Marchiolli and D. Galetti for lengthy and 
fruitful discussions. 



Appendix A. Approximate unitarity 

In this appendix we may define the conditions in which we can consider an operator 
of the following form 



as an "approximately" unitary operator. It's hermitian conjugate is directly calculated 



iV,-l 



N-l 




(A.1) 



JV|-1 N-l 



(0= E M)M+ E l»;i>(»;i 
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Ni-l Nt-1 N t -1 N-1 

z u ml(0= E E |n;0(n;iKi)(m ; ^|+E E l«;0(n;iKi>(m;i| 

n— m—0 n— ra—Ni 

N-1 Nt-1 N-1 N-1 

E E |n;l)(n;l|m;l)(m;C|+ E E |n; l)<n; l|m; l)(m; 1| 

n—Ni m—0 n—Ni m—Ni 

In order to study the situations in which the above operator is equal to the identity 
operator, the first needed criterion is Vf and are "approximately" orthogonal. We 
establish this condition as 

|(m;£|n;£}| 2 < 1(T 4 , for n<Ni<m. 

We now have 

N,-lNi-l N-1 N-1 

S„(04(0=E E \n;Z)(n;l\m;l){m;Z\ + E E M)(n; l|m; l)(m; 1|. 

n— m—0 n—Ni m—Ni 

The second term is a different representation the identity operator of [15], whereas 
the first one is an "approximate" representation of the identity operator of V ; as long 
as (n; l\m; 1) ~ (n; £|m; £). To establish objectively this criterion we define that when 

|(n;lKl)| 2 -|(n;CKO| 2 <10- 4 

the overlap (n;£|ra;£) is assumed not to depend on £, so under this approximations 

= h + h = i. 

The conditions emerging from the relation Sjj(£)S„(£) = 1, combined with the ones 
above can finally be summarized in the conditions 

\(m;£_\n;£_)\ 2 < 1CT 4 , for n < N t < m, 

|(n;l|m; 1)| 2 - \(n^\m;0\ 2 < 10" 4 , 

as the relations involving the dual vectors are consequences of those. 
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